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ABSTRACT 

In 2008 AGILE and Fermi detected gamma-ray flaring activity from the 
unidentified EGRET source 3EG J1236+0457, recently associated with a flat 
spectrum radio quasar (GB6 J1239+0443) at z=1.762. The optical counter- 
part of the gamma-ray source underwent a flux enhancement of a factor 15-30 
in 6 years, and of ~10 in six months. We interpret this flare-up in terms of 
a transition from an accretion-disk dominated emission to a synchrotron-jet 
dominated one. 

We analysed a Sloan Digital Sky Survey (SDSS) archival optical spectrum 
taken during a period of low radio and optical activity of the source. We esti- 
mated the mass of the central black hole using the width of the C iv emission 
line. In our work, we have also investigated SDSS archival optical photomet- 
ric data and UV GALEX observations to estimate the thermal-disk emission 
contribution of GB6 J1239+0443. 

Our analysis of the gamma-ray data taken during the flaring episodes indi- 
cates a flat gamma-ray spectrum, with an extension of up to 15 GeV, with no 
statistically-relevant sign of absorption from the broad line region, suggest- 
ing that the blazar-zone is located beyond the broad line region. This result 
is confirmed by the modeling of the broad-band spectral energy distribution 
(well constrained by the available multiwavelength data) of the flaring activity 
periods and by the accretion disk luminosity and black hole mass estimated 
by us using archival data. 

Key words: galaxies: active - galaxies: quasars: general - galaxies: individual: 
GB6 J1239+0443 - radiation mechanism: non thermal 



1 INTRODUCTION 

Blazars are a sub-class of active galactic nuclei, emitting from radio to TeV energies. They 
are subdivided in two main categories: Flat Spectrum Radio Quasars (FSRQ) and BL Lac- 
ertae (BL Lac) objects. FSRQs are characterized by a flat radio spectrum in the GHz range 
(with spectral index a ^ 0.5, where the flux density is S v oc v~ a \ and strong and broad 
emission lines (with equivalent width >5 A). BL Lac objects, on the other hand, have no or 
weak emission lines with equivalent width <5 A. 

Blazars continuum emission originates from a relativistic jet aligned wi th the line of sight. 



Their Spectral Energy Distribution (SED) shows a double humped shape (jUrry and Padovani 



19951 ). with a low energy peak lying between IR and X-rays, and an high energy peak in the 
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MeV-TeV band. 

The low energy region of blazar spectra is associated with the synchrotron emission coming 
from the jet relativistic electrons. The high energy region can be modeled through inverse 
Compton emission (leptonic models), with seed photons coming from an external region 
(e.g. the accretion disk, the dusty torus), eventually reprocessed by the broad line region, 
or the hot corona, or from the synchrotron process itself (synch rotron self-Compton, o r 



SSC). A detailed descrip ti on of leptonic mod els can be found in 



Marsher fc Bloom 



(119921 ): 



Sikora et al. 



(1199 



Maraschi et al. 



(|1992h : 



T he high energy region can also be mo deled with hadronic scenarios (jMiicke and Protheroe 



2001 



iMiicke et al. 



2003 



Bottcher 



20071 ). where the very high energy protons of the jet are 



radiatively important. The accelerated protons produce gamma ray emission through pro- 
ton synchrotron emission, the decay of neutral pions, and synchrotron emission produced 
by secondaries. 

The location of the so-called "blazar- zone" , i.e., the spatial location of t he blazar SED - peaks 
and gamma-ray emitting region, in blazars is still a matter of debate. 



Sikora et al. 



d2008|) 



proposed that the blazar-zone is located at 3-9 pc from the cen tral engine for the outbu rst of 



Ghisellini et al. 



( 120071 ) indi- 



3C 454.3 (a bright FSRQ) occurred in 2005. For the same flare, 
cated, instead, a dissipation region at 0.5-0.8 pc from the central black hole (BH). From the 
combined study of time-dependent polarimetric radio image s at 43 and 86 GHz, t he optical 



polarimetry, and radio, optical, X-, gamma-ray light curves, 



Jorstad et al. 



( 120101 ) proposed 



that the low and high energy emission is located near the 43 GHz core, at a distance of the 



order of tens of parsecs from the BH for 3C 454.3. A simi 



BL Lac objects OJ2 8 7 and AO 0235+164 (lAeudo et al. 



Tavecchio fc Mazin 



ar investigation, performed for the 
20 11 al ibi), led to similar results. 



(120091 ) established that gamma-rays emitted inside the broad line re- 
gion (BLR) are absorbed at E> 20 GeV/(l + z) due t o the 77 interaction with the BLR 



photons (internal absorption). 



Poutanen fc Stern 



(120101 ) refined this result, and claimed in- 



ternal absorption features at E> 5 GeV/ (1+z) and at E> 20 GeV/ (1 + z) in the gamma-ray 
spectrum of 3C 279, 3C 454.3 , PKS 1510-08, and a few oth er FSRQ. 



Within the leptonic scenario, 



Ghisellini fc Tavecchio! ( 120091 ) show that the contribution of 



external photon fields, including contributions from the BLR and dusty torus to the inverse 
Compton emission can be parametrized as a function of the accretion disk luminosity an d 
the dissipation distanc e of the emitting blob from the BH (jGhisellini fc Tavecchio! 120091 ) . 



Ghisellini et al. 



(120 lOl ) modeled the SED of the gamma-ray loudest blazars as being emit- 
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ted at 0.01-0.5 pc from the BH. 



Using multiwavelength observation s of the blazar GB6 J1239+0443, we will apply the 



models of 



Ghisellini fc Tavecchiol ( 120091 ) to further investigate the location of the blazar-zone. 



GB6 J1239+0443 was an uniden tified gamma-ray sour c e of the Virgo region (3EG J 12 36+045 7) 



detected with low significance (IHartman et al. 



1999 



Casandjian fc Grenier Il2008|) by the 



EGR ET gamma-ray telescope (operating in the 30 MeV - 30 GeV energy range, see 



Esposito et al. 



19991 ). In recent years the gamma-ray source has shown two episodes of remarkable high en- 



ergy activity: at the b eginning and at the e nd of 2008, when it was detected by the AGILE 



respectively; t 


len n 


(Nolan et al. 


2012) 



hen named AGL J1238+0406 ( Pittori et al. 



20091), and 2FGL J1239. 5+0443 



the association of the unidentified gamma-ray source with BZQ 1239+0443 a flat spec- 
trum radio quasar (FSRQ) in cluded into the secon d edition of the Roma-BZCAT Multi- 
Frequency Catalog of Blazars (iMassaro et al. II2010I ). The optical counterpart of this source 
is SDSS J123932. 75+044305. 3, located at z=1.762, and the radio counterpart is named 
GB6 J1239+0443. In the following sections, we refer to this object as GB6 J1239+0443, 
its radio source name. 

Here we present the results of an analysis of multifrequency data simultaneous to the AGILE 
campaign on the Virgo field, and to the follow-up carried out after the Fermi-LAT detec- 
tion and localization. By analysing the archival data, we estimate some fundamental physical 
properties of the source such as the accretion disk luminosity and the supermassive black 
hole (SMBH) mass. In this way we can obtain a consistent picture of the source emission 
in the framework of leptonic models of blazars for periods of both low and high emission 
activity. In Sections 2 and 3 we will describe the multi-wavelength campaigns related to 
this source. In Section 4 we will report on the archival data. In section 5 we will present 
our results consisting in the determination of the BH mass, the gamma-ray light curve and 
spectrum, and the SED modeling. 
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Table 1. Observing campaigns of the Virgo field during high gamma-ray activity periods from GB6 J1239+0443. 



campaign observatory 



observing period 



from 2007 Dec. 16 17:14 to 2007 Dec 23 02:18 UT 

AGILE from 2007 Dec. 24 07:12 to 2007 Dec. 30 23:03 UT 

from 2008 Jan. 04 13:35 to 2008 Jan. 08 11:06 UT 
from 2007 Dec. 19 18:08 to 2007 Dec. 22 06:44 UT (revolution 633) 
INTEGRAL from 2007 Dec. 25 17:39 to 2007 Dec. 28 06:27 UT (revolution 635) 

from 2007 Dec. 31 17:13 to 2008 Jan. 03 04:00 UT (revolution 637) 



Fermi 



Dec. 2008 - Jan. 2009 



Swift 



from 2009 Jan. 02 17:47 to 2009 Jan. 04 20:21 UT* (total exposure 4.7ks) 



KANATA 



2009 Jan. 02 18:14 UT for 1.5ks 



2009 Jan. 03 19:41 UT for 1.5ks 



Two pointings on 2009 Jan. 02 (exposure 2562 s and 380 s), and 5 pointings on 2009 Jan. 04 (exposure 597 s, 537s 

577 s, 537 s) 



537 s, 



2 GAMMA-RAY OBSERVATIONS AND RELATED MULTIFREQUENCY 
CAMPAIGNS 

Th e A GILE-GRID gam ma-ray telescope (operating in the 30 MeV - 50 GeV energy range, 



sec 



Tavani et al. 1 120091 ) performed two observing campaigns of the Virgo field. The first 



campaign included 3 observations from 2007 December 16 to 2008 January 8. There were 
3 simultaneous observations (revolutions 633, 635, 637) with the wide-field instruments 
aboard the INTEGR AL mission (operating in optical, hard X-, and soft gamma-rays, see 



Winkler et al. 



20031). During this camp aign AGILE detected high gamma-ray activity from 



GB6 J1239+0443 (jPacciani et al. II2009I ). In the following sections, we will refer to this time 
period as period A, and to the related multiwavelength campaign as campaign A. 
AGILE also observed the Virgo field from 2009 June 4 to 2009 June 15, but there were 
no simultaneous observations with wide-field instruments operating at other wavelengths. 
GB6 J1239+0443 was undetected during this observation. 



The Fermi-LAT gamma-ray telescope (20 MeV - 300 GeV, see 



Atwood et al. 



20091 ). op- 



erating in all-sky survey m ode since 2008 August 4, detected high gamma-ray activity at 



the end of December 2008 (ITramacere et al 



200 9J), a nd triggered optical-UV/X-ray obser- 



vations with the Swift satellite ( jGehrels et al. 1120041 ). starting from 2009 January 2 for a 
total exposure of 4.7 ks. Optical V band photometry and polarization measurements were 
al so made on-ground with the KANATA telescope on 2009 January 2 and on 2009 January 
3 (jlkejiri et al. II2009I ). In the following sections, we will refer to the observations collected 
during this period as period B, and to the related multiwavelength campaign as campaign 
B. A summary of the observations is provided in Table [TJ 
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3 DATA ANALYSIS 



3.1 AGILE-GRIT) data 



Level-1 AGILE-GRID data for campaign A were analysed using the AGILE Standard Anal- 



ysis Pipe 



ine (BUILD20) and the AGILE Scientific Analysis Package, based on the likelihood 



method (IMattox et al. 



19961 ). Albedo photons were rejected by applying a cut at 85° cen- 



tered on the Earth. We selected well-reconstructed g amma-rays by applyin g the FM3.119 



filter, calibrated in the 100 MeV - 3 GeV energy band ( ICattaneo et al. 



2011). All the events 



collected during the passage in the South- Atlantic Anomaly were rejected. Counts, exposure 
and Galactic background gamma-ray maps were created with a bin-size of 0°.l x0 o .l for E > 



100 M eV. We detected a source (A GL J1238+0406 in the AGILE catalog, see 



Pittori et al. 



2009, and 



Verrecchia et al. 



rameter, see 



Mattox et al. 



2011) with a significance of ~6 (as measured by the yTS pa- 
1996^ , located at a2ooo=190.25, <52ooo=4.40, with an error radius 



of (33 + 6) arcmin (statistical error at 95% confidence level -C.L.-, and systematic error, 
respectively), by integrating the GRID data for 4 days between 2008 January 4 13:35 and 
2008 January 8 11:16 (within the campaign A). The source was positionally consistent with 
GB6 J1239+0443. A multisource maximum likelihood analysis was performed to extract the 
source flux and position taking into account nearby sources 3C 273, 3C 279, and 4C 04.42 
(for which we obtained yTS > 1 from a preliminary analysis of the observing campaign). 
For AGL J1238+0406 we obtained a flux of (62±9)xl0~ 8 ph cm~ 2 s" 1 (E > 100 MeV) and a 
photon index 1.92±0.14. The integration of the first week of observations with AGILE gave 
no detection at the position of GB6 J1239+0443, resulting in an upper limit of 21xl0 -8 ph 



— 2 —1 

cm s . 



3.2 Fermi-LAT data 

We analysed the Fermi Large Area Telescope (Fermi-LAT) data for campaign B with the 
standard Fermi Science tools v9r23pl, following the prescriptions in the online documenta- 
tion 0. We used the Pass 7 response functions (P7_V6). In particular, we selected events of 
event class 2, suitable for point-like sources, and we filtered out photons from the Earth's 
limb with a cut at 100° in the zenith angle. We performed the unbinned likelihood analysis 
inside a region of radius 15° around GB6 J1239+0443 to derive the source flux. We took 



http: / / fermi . gsfc. nasa.gov / ssc / data / analysis / documentation / 
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into account the diffuse backgrounds, which were modeled using gal_2yearp7v6_v0 for the 
Galactic diffuse emission and iso_p7v6source for the extragalactic isotropic emis sion models 



and a ll the 51 point-like gamma-ray sources in the second Fermi-LAT catalog ( iNolan et al. 



20121 ) within a slightly larger radius of 20° from GB6 J1239+0443 (we considered a larger 



radius due to the PSF width). For each source we used the model specified in the second 
Fermi-LAT catalog. For the sources within 10° from GB6 J1239+0443, we kept free all the 
spectral parameters in the fit. For the sources within an annulus of internal radius 10° and 
external radius 15° we kept free only the parameters related to the flux normalization, and 
all the other parameters were fixed to the values reported in the second Fermi-LAT catalog. 
For the sources outside 15° from GB6 J1239+0443, we fixed all the spectral parameters to 
the values reported in the catalog. This is a standard procedure for the analysis of Fermi- 
LAT data, implemented with the python routine make2FGLxml.py3 (contributed software 
by T. Johnson) 

We proceeded with the analysis for energies only above 300 MeV, in order to process data 
with a smaller point spread function and reduce background gamma-rays from 3C 273, a 
bright and soft gamma-ray source (r ~ 2.6) located at ~4° from GB6 J1239+0443, since 
the 68% (95%) containment radius for Fermi-LAT at normal incidence is 4°. 5 (10°) at 100 
MeV @. Our study showed that this choice had a negligible effect on the signal significance 
of GB6 J1239+0443. 

We used the gtfindsrc tool to locate the gamma-ray source. By integrating data for 1 month 
centered around the peak flux (2008 December 29 16:00 UT) in the band 300 MeV - 20 GeV, 
and by using photons converted both in the front- and back-section of the Fermi-LAT , we 
obtained a detectio n of a source (2FGL J1239. 5+0443 in the second Fermi-LAT catalog, see 



Nolan et al. 



20121 ) with \/TS=20 located at a 2 ooo=189.897, 5 2 ooo=4.718 and an error radius 



of 9 arcmin. The source was positionally consistent with GB6 J1239+0443. We obtained a 
gamma-ray photon index of 2.15 ± 0.11. 

From the analysis in the 300 MeV - 20 GeV range, we also obtained a detection with 
VTS ~18, a flux of (23±3)xl0^ 8 ph cm~ 2 s -1 , and a photon index of 2.21 ± 0.15, when 
keeping the integration time within only 4 days centered at 2008 December 29 16:00 UT 
(the peak flux). During this integration period, the source was detected up to the energy 



2 
3 
4 



both available on the URL http://fermi.gsfc.nasa.gov/ssc/data/access/lat/BackgroundModels.html 

available on the URL http://fermi.gsfc.nasa.gov/ssc/data/analysis/user/ 

http: / /www. slac.stanford.edu/exp/glast/groups/canda/archive/lat_Performance. htm 
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interval 10-20 GeV, for which we obtained a detection with y/TS = 5.8. To compute the 
upper limits needed to build the source light curve and spectra, we used the UpperLimits 
python function provided with the Fermi Science tools. 



3.3 INTEGRAL/IBIS data 

The INTEGRAL/IBIS (operating in the 17 - 400 keV energy range, see 



Ubertini et al. 



20031 ) data for campaign A were processed using the OSA software version 8.0. We searched 



for the source starting from the images accumulated in the 20-40 keV band for revolutions 
633, 635, 637 (simultaneous with A GILE observation of the Virgo field, campaign A). IBIS 
did not detect the source. We derived a 3 sigma upper limit of 1.7 mCrab for each revolution 
(200 ks exposure). 



3.4 Swift-XRT data 

The Swift-XKT (X-ray Telescope, operating in the 0.2 - 10 keV range, see 



Burrows et al. 



20051 ) data for campaign B were processed using the most recent available calibration files. 



We made use of Swift Software version 3.5, FTOOLS version 6.8, and XSPEC version 12.5. 
The observations were obtained in photon counting mode, with a total integration of 4.7 ks. 
The mean source count rate was (2.58±0.23)xl0~ 2 cps. We extracted the spectrum using 
a photon binning ratio that ensured more than 20 photon counts per energy bin. We fitted 
the X-ray data with an absorbed power la w, fixing the absorption to the galactic value of 
1.85xl0 20 cm" 2 ( IDickev fc Lockman Hl99(I ). We obtained a photon index of 1.42±0.25 (90% 
C.L.). The estimated flux in the range 2-10 keV was (8.8±2.7) x 10~ 13 erg cm" 2 s" 1 (68% 
C.L.). 



3.5 Swift-UVOT data 



Swift-XJVOT (Ultra- Violet /Optical Telescope, see lRoming et al. 1120051 ) data from each ob- 
servation sequence of period B were processed by the standard UVOT tool uvotsource 
using the same version of the Swift software as for the XRT analysis. An extraction re- 
gion of radius 5 arcsec centered on the source and a background region of radius 13 arcsec 
located at a 20 oo =12 h 39 m 29!66, 5 20 oo=+04 42'34" 2 (at least 27 arcsec far away from any ob- 
ject in the NED database) were used. Magnitudes are expressed in the UVOT photometric 
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system (IPoole et al. 1120081 ). We obtained m[/=16.27±0.03 for GB6 J 1239+0443 (extin ction 
corrected using the mean Galactic interstellar extinction curve from iFitzpatrick Ill999l ). 



3.6 Kanata optical data 

The optical photometry was performed using TRIS PEC (a simultaneous o ptical and near- 
infrared imager, spectrograph, and polarimeter, see Iwatanabe et ah" 2005 ) attached to the 
Kanata 1.5-m telescope at Higashi-Hiroshima Observatory on 2009 January 2 at 18:14 and 
2009 January 3 at 19:41 UT (period B). The observations were performed in polarimetry- 
mode with a narrow aperture mask of 1.5' width. The total exposure was 1476 s per night. 
The observations were pipeline-reduced, including bias removal and flat-field corrections. 
We derived the V-band magnitude from differential photometry with a nearby reference 
star at a 2 ooo =12 /l 39 m 30 ! !ll, 5 2 ooo=+04°39'52 , . / 6 of which the V magnitude of 14.095 was 
deduced from the q' and r' data in the Q th release of the Sloan Digital Sky Survey (SDSS, 



Adelman-McCharty et al. 



20081 ). 



3.7 INTEGRAL/ OMC data 

We analysed the OMC (Optical Monitoring Camera, equipped with a Johnson V filter, see 



Mas-Hesse et al. 



20031 1 data collected during period A with the OSA software version 8.0. 



Due to the dithering mode of INTEGRAL observations, the field around the source was 
observed with OMC in the science windows 2, 50, 52, 78 of revolution 633, in the science 
windows 49, 51, 77 of revolution 635, and in the science windows 49 and 51 of revolution 
637. To minimize readout noise contribution to the signal to noise ratio (S/NR), we selected 
snapshots of 200 s. Typically 6-10 snapshots of at least 200 s integration were recorded 
for each science window. Due to telemetry limits, the whole frame of the OMC field of 
view cannot be telemetred, but only sub-frames around sources from a reference catalog. 
GB6 J 1239+0443 is included in one of the predefined sub-frames. We detected the source 
with S/NR=4 and my=17.5±0.3 with an aperture of 3x3 pixels, averaging the 3 revolutions. 
The investigation of previous INTEGRAL observations of the field resulted only in upper 
limits (my > 18.1-19.1 depending on the exposure of the observations). The Mosaic image 
from the individual snapshots taken during the multiwavelength campaign (obtained using 
the iraf package) is shown in Fig. [IJ The pixel size of the mosaic image was 1/3 of the 
OMC pixel to reduce spurious effects in the images shift. A cosmic ray signal was filtered 



Flaring and low activity periods of GB6 J1239+0443 11 




Figure 1. Mosaic image obtained with OMC simultaneous to the AGILE detection of gamma-ray flare of GB6 J1239+0443. 
The pixel size of the mosaic image is 1/3 of the original OMC pixel. The label AGN1 is for GB6 J1239+0443 
(SDSS J123932. 75+044305. 3), AGN2 for SDSS J123936. 52+044249.3, AGN3 for SDSS J123945. 06+044431.4. The circle 
marked C.R. indicates the filtered region where a cosmic-ray event was collected during revolution 635. 

out (affecting the combined image around a2ooo=12 39 32.5, <52ooo=+04 44 44.5), for science 
window 77 of revolution 635 only. This region is marked as C.R. in Fig. [TJ 



3.8 VLA archival data 

We retrieved archival Very Large Array (VLA^I) D-array 43 GHz data of GB6 J1239+0443, 
observed on 2001 November 5 for 65 s. The VLA data were reduced and analysed using 
standard routines implemented in the Astronomical Image Processing System (AIPS) pack- 
age. The sources 3C286 (1.45 Jy) and 1254+116 (0.43 Jy) were used as flux and phase 
calibrators, respectively. 

An almost unresolved source has been detected at 43 GHz with position aj2ooo=12 h 39 m 32!75±O f 501 
and 5j2ooo = 04°43 / 05'. / 2±0'.'l, a peak flux density of ~71.4±3.6 mJy, and an integrated flux 
density of ~81.2±4.1 mJy. 



4 THE ARCHIVAL DATA 

In order to model accurately the SED (sections 5.2 and 5.4), we used archival data of 
GB6 J1239+0443 from the SDSS. Photometry with the filters u, g, r, i, z was performed on 
2001 March 15. In particular the SDSS archive reports m u =20.62±0.06 and m g =20.47±0.03 
(corresponding to mu ~19.9, and my ~20.5). An optical spectrum (370-920 nm) was ob- 
tained on 2002 May 13. 



5 The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under cooperative 
agreement by Associated Universities, Inc. 
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Table 2. Summary of archival optical and radio observations of GB6 J 1239+0443. 



date 



proj ect / observatory 



measurement description 



2001 March 15 
2001 November 5 
2002 May 11 
2002 May 13 
2007 January 18 
2007 April 17 & May 13 
2009 January 30 

2009 December 10 

2010 January 3-9 



SDSS 
VLA 

Metsahovi 
SDSS 

UKIDSS large Area Survey 

GALEX 

VLBA 

VLBA 

Planck 



optical photometry 
43 GHz 
22 GHz 
Optical Spectrum 
near infrared photometry 
near UV photometry 
15 GHz 
15 GHz 
30 - 217 GHz 



The near infrared photometry of the source was performed by using the data of the UKIDSS 
(UKIRT Infrared Deep Sky Survey) LAS (Large A rea Survey) of 2007 January 18. 



GALEX (an orbiting ultraviolet space telescope, see lMartin et al. 1120051 ) observed the source 
with NUV and FUV filters (with bandpass centered at 230 and 150 nm respectively) between 
2007 April 17 and 2007 May 13. 

The source was observed by VLA at 43 GHz on 2001 November 5, resulting in a low flux 
(we reported the analysis of this observation in the previous section). 

A weak detection at 22 GHz was obtained by the Metsahovi observatory on 2002 May 11 
( ITerasranta et al. IboOJ ), they reported a flux of 0.22±0.04 Jy. 

Following the Fermi-LAT detection of a gamma-ray flare, the source has been added to the 
MO J AVE 3 sample, and was observed twice by VLBA in 2009 at 2 cm: on 2009 January 30 
(one month after the gamma-ray flare) and on 2009 December 10. The source was found in 
the Planck legacy archive v0.5, with detections at 30, 100, 143, 217 GHz with observations 
on 2010 January 3, 2010 January 6, 2010 January 9, 2010 January 7, respectively. We also 
added the 147 GHz data for which the detection is flagged as extended. At last we included 
in the SED the detection in the ROSAT all Sky Survey. A summary of the archival obser- 
vations is given in table [2l while the archival radio and optical SED is plotted in figure [21 



6 MOJAVE (Monitoring Of Jets in Active galactic nuclei with VLBA Experiments) is a long-term program to monitor radio 
brightness and polarization in jets of AGN. 
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Figure 2. Top panel: archival radio observations of GB6 J1239+0443. Yellow box data was taken by VLA in 2001, Green data 
was taken by Metsahovi in 2002. Cyan diamond taken by VLBA in December 2009. Cyan stars are for Planck taken in 2010. 
Purple data represent the higher and lower fluxes measured at 8.4 GHz and reported in literature (from Mojave database). 
Bottom panel: archival optical observations of GB6 J1239+0443, all data are Galactic absorption corrected. Yellow squares are 
SDSS photometry taken in 2001; the green line is the SDSS optical spectrum taken in 2002 (squares of the same color with 
superposed black stars represent the simulated SDSS photometry evaluated from the observed spectrum). Light green squares 
show near-IR photometry taken by UKIDSS-LAS at the beginning of 2009, and blue squares are the GALEX data taken in 
mid-2007. The dashed and dot-dashed lines represent two disk-emission models; see text for details. 
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5 RESULTS 



5.1 Disk luminosity and black hole mass determination from the archival 
optical/UV photometry 

In spite of the large amount of data available for GB6 J1239+0443 the source has not been 
studied in detail before. 

The optical SPSS spectrum was ta ken simultaneously with the 22 GHz Metsahovi radio 
observation (ITerasranta et al. II2005I ). We note that the 22 GHz flux is a factor 3-4 lower 
than the Planck data of 2010 (light blue data in top panel of fig. [2]). From the optical spectrum 
shape (bottom panel of fig. [2]) and from the faint radio emission, a reasonable assumption is 
that the optical emission in M ay 2002 is accretio n-disk-dominated. We evaluated the non- 
thermal dominance (NTD, see 



Shaw et al. 



2012[ ) of the optical spectrum making use of the 



C IV line luminosity, and the continuum luminosity extrapolated at 1350A rest-frame. 

We obtained NTD=l.l±il where the ma j or con tribution to the errors comes from the 
uncertainties in equation (3) of 



Shaw et al. 



fl2012f ). and NTD=1.0±%i making use of the 



Mg II line luminosity and the continuum luminosity at 3000 A rest-frame, indicating that 
the optical spectrum taken in 2002 is typical of a thermal emission dominated source. We 
also observe that the optical photometry taken by SDSS in March 2001 results in a similar 
shape and lower flux (~ 30% lower), and that the radio observation by VLA in November 
2001 shows the source in a rather radio faint state. Therefore we also make the assumption 
of accretion disk dominated emission for the optical photometry taken in March 2001. 



We compare the optical data taken by UKIDSS-LAS in January 2007, the SDSS spectrum 
taken in May 2002, and the optical photometry taken in March 2001. These observations 
taken at different epochs seem to demonstrate three different emission states: a jet-dominated 
emission state (UKIDSS-LAS), an accretion-disk-dominated state (SDSS photometry), and 
an intermediate state (SDSS optical spectrum) that seems to be dominated by accretion-disk 
emission with a possible minor contribution from the jet. We also note that the GALEX 
observations give a factor of ~10 lower flux with respect to the Swift/UVOT observation 
taken during the 2009 January flare. We cannot establish whether the UV emission detected 
by GALEX is due to the disk alone or whether there is also the contribution of the high en- 
ergy tail of synchrotron emission. In the latter, unfavoured, case, the UV data from GALEX 
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result in at least an upper limit on the disk emission. 



With these considerations, and assuming the disk emission is not variable on year 
timescales, we attem pted to model the SPSS an d GALEX photometry with a Shakura Sun- 
yaev accretion disk (IShakura fc Sunyaev Ill973l ) around a non-rotating black h ole (dashed 



curve in the bottom panel of Fig. [2]) . We modeled the disk emission as proposed by 
(120091 ). with inner radius of 3 Schwarzschild radii (Rs), and outer radius of 500 Rs- We fit 
the model to the SDSS+GALEX data keeping free the parameters of disk luminosity (L^) 
and Rs- We obtained a disk luminosity of ~8.9xl0 45 erg s _1 , and i?s~2.4x 10 14 cm, corre- 
sponding to a maximum emitting temperature of ~5.4xl0 4 K, a BH mass of ~8xl0 8 solar 
masses, an accretion rate of ~9 M ?/ _1 (assuming an accretion efficiency e accr =0.1), and 
then an Eddington ratio Re<m = L d /L Edd ~ 50%. 



Ghisellini fc Tavecchio 



We note that one observation by SDSS and one observation by GALEX show a flux 
higher than the one expected according to our thermal emission model, but the filters used 
include quasar line emission from Mgll (the SDSS i filter) and Lya (the GALEX NUV 
filter); see Figure El We also fit a disk model assuming the near-UV (GALEX) photometry 
is dominated by jet emission. In this case, the model assumes an accretion disk with only 
a minor contribution from jet emission and is fit to the SDSS data alone. In this scenario 
(dot-dashed curve in the bottom panel of Fig. [2]), a model with a disk luminosity >5.4xl0 45 
erg s _1 is required. We are aware that our model fits are not unique and that there are no 
simultaneous radio observations to strongly validate one fit over the other, but our preferred 
model reproduces the ob served fluxes and is con sistent with the scenario usually proposed for 

where hard optical/UV spectra are usually 



Ghisellini et al. 



2011 



FSRQs (see for example 
modeled as accretion disk dominated emission). 



5.2 Black Hole Mass Determination from Archival Optical Spectrum 



The BH mass of GB6 J1239+0443 could be es timated better by using t h e sing le-epoch BH 



mass scaling relationship for C IV derived from iVestergaard fc Peterson I (120061 ) and applied 
on the archival optical spectrum for this FSRQ. Unfortunately, the SDSS spectrum has 
a rather low S/NR (~ 3 per pixel in the continuum near the CivA1549 emission line). 
Penney et al. I (120091 ) show that line widths and thus single-epoch BH masses measured 
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Figure 3. The C IV line in the archive optical spectrum of GB6 J 1239+0443. The red continuum line is the fit with a sixth-order 
Gauss-Hermite polynomial. 

from low S/NR data have relatively larger systematic uncertainties than those measured 
from high quality data. I n addition, unrecognized absorption is a particu l ar concern for 



low q uality C IV data (see IVestergaard fc Peterson 



2006 



Assef et al. 



2011 



Penney et al. 



201ll ). Nonetheless, this is the only optical spectrum during a low state currently available. 
Therefore we use this spectrum to measure the line width and determine the mass. We em- 
ployed two line width measurement methods and quote conservative uncertainties on the line 
width taking into account the low quality of the data. Our methods for measu ring the C IV 



Assef et al. 



l ine w idth and uncertainties closely follow the prescription "A" described by 
(120111 ). and we therefore refer the reader to this work for details. After subtracting the lin- 
early fit continuum, based on the windows shown in Fig. [3j we measured the FWHM of C IV 
both directly from the data (FWHM = 2860 ± 910 kms -1 ) and from a sixth-order Gauss- 
Hermite polyn o mial fi t to the line profile, as shown in Fig. [3] (FWHM = 4710 ± 390 kms -1 ). 
Penney et al. I (120091 ) show that direct measurement of the FWHM from low S/NR data 
systematically underestimates the line width, while measurement from a fit referring to the 
same data can overestimate the same width. We adopt a conservative approach and take 
the mean of these two width measurements and assume the quadrature sum of the uncer- 
tainties. The adopted FWHM measurement becomes FWHMciv = 3800 ± 1000 km s -1 . 
We then measure the mean continuum luminosity in the continuum window near rest-frame 
1450A to be AL1450 = (3.47 ± 0.44) x 10 45 erg s _1 , after correcting for Galac t ic ext inction. 



We then evaluate the mass using equation (7) from IVestergaard fc Peterson I (120061 ). which 



is also Equation (6) from 



Assef et al. 
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t is worth noting that the SPS S spectrum 



(2011). 



does not extend to rest-frame 1350A; however, 



Vestergaard fc Peterson 



(120061 ) argue that 

the 1450A luminosity can be substituted without penalty for the 1350A luminosity, as we 
have done h e re. W e estimate the BH mass of GB6 J1239+0443 to be 4.3l|;| x 1O 8 M . 



Assef et al. 



(1201 if ) find a correlation between the ratio of the C iv-to-Balmer mass estimates 
and the UV-to-optical luminosity ratio. Since this correlation is based on the ratio of the 
mass estimates, barring further investigation into the source of this correlation, it is unclear 
whether it is the Civ-based or Balmer-based mass estimates, or both, to be the source of 
the bias. Regardless of origin, Assef et al. found that, when this correlation is removed, 
and when they arbitrarily choose to correct the Civ-based masses, the corrected masses 
are highly consistent with the measured Balmer-line-based mass estimates (the scatter in 



the corrected Civ vs. Balmer mass estimates is reduced 
to the uncorrected mass estimates; see 



Assef et al. 



Dy a factor > 2 when compared 



20111 ). Balmer-based mass estimates 



are generally more accepted in the literature because they are relativel y better calibrated 



with direct mass measuremen t s from reverberation mapping (see, e.g 



Collin et al. 



200 



Vestergaard fc Peterson 



2006 



Penney et al. 



* 



20091 ). At this point, however, it is impossible 



to state which of the two mass measures is actually more accurate. In particular, we must 
consider that host galaxy starlight can significantly contaminate the optical luminosity with 
which Balmer-based masses are estimated, yet there may be evidence of non-virial motions 



Richards et al. 



20021 ). For GB6 J1239+0443, we fit a powerlaw con- 



from Civ (see, e.g. 

tinuum to the full wavelength extent of the SPSS spectrum and extrapolate to rest-frame 
5100A to estimate the res t-frame optical lum inosity to be AL 5100 = (2.64 ± 0.33) x 10 45 erg 



s _1 . Using equation (8) of lAssef et al. I (1201 ll ). with the coefficients based on their prescrip- 
tion A, we can then calculate a corrected Civ-based BH mass of 5.3i 3 ; 3 x 10 8 M Q . 
The two mass estimates (based on the C IV broad line width and on the thermal continuum) 
are in qualitative agreement. 



5.3 Gamma-ray light curve 

We created a gamma-ray light curve for the source from both the AGILE pointing and the 
Fermi survey, as shown in Fig. HI The AGILE data were integrated with bin sizes of about 
6.5, 6.5, and 4 days, due to gaps in the observation (see Tabled]). The Fermi-LAT data were 
integrated with binsizes of 1 and 7 days. The flux reported is in the 100 MeV - 3 GeV range 
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Figure 4. The gamma-ray light curve for GB6 J1239+0443 obtained with AGILE during campaign A (diamond symbols and 
right vertical scale for the gamma-ray flux) and Fermi (left vertical scale for the gamma-ray flux). For the Fermi-LAT data we 
report also the photon index evaluated in the 300 MeV - 20 GeV energy band. The upper panel reports Fermi— LAT data with 
a typical binsize of 1 day. The lower panel reports Fermi— LAT data with 7 days integration. AGILE data were integrated with 
bin sizes of about 6.5, 6.5, 4 days, due to gaps in the observation (see Table [TJ- The orange bands represent the INTEGRAL 
campaign (campaign A); width is in scale. The green bands represent the Swift campaign (campaign B); width not in scale. 



for AGILE. As discussed in section 3, the Fermi-LAT data were analysed between 300 MeV 
and 20 GeV in order to reduce the contamination from the nearby 3C 273 at lower energies. 
In Fig. [H we also report the gamma-ray photon index in the 300 MeV - 20 GeV as obtained 
by the Fermi-LAT data assuming a power-law spectrum. 



b ev aluate the flare duration, we made use of the procedure described by 



Abdo et al. 



((20101). From the light curve with time bin of 1 day, we obtaine d a duration (defined as 



"^ T}al1 ) of 7 days and an asymmetry -0.3 (see lAbdo et al. 



20101 for details, T rise and Tf a u 



are the rising and falling time respectively). The light-curve with time bin of 7 days shows 
more than a single relative maximum; therefore, the fit with a simple curve is not feasible. A 
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Figure 5. Gamma-ray spectrum of GB6 J1239+0443 integrated for one month around the peak flux (black data). The spectrum 
integrated for 4 days around the flare is reported too, multiplied by a factor 2 (cyan data). The fitted models discussed in the 
paper are reported. The continuous lines are for the models, the binned curve are for models weighted over each energy bin. 
The power law model is reported in orange, the power law model + gamma-gamma absorption is reported in green. 



rough definition of the gamma-ray activity period could instead be the timespan for which 
Fermi-LAT detects gamma-ray emission from the source. Assuming temporal bins of 7, 15, 
and 30 days, we found that Fermi-LAT detected gamma-rays from the source for at least 
11 weeks. 



5.4 The Fermi- -LAT gamma-ray spectrum for the high activity period of 
campaign B 

The gamma-ray spectrum obtained by Fermi-LAT data, integrated for one month cen- 
tered around the peak flux (2008 December 29 16:00 UT), is reported in Figure |5] together 
with the spectrum integrated for four days. The 30-day integrated spectrum has been fit 
first with a power-law model. Because absorption is expected for a blazar-zone originating 
near the central source (e.g., within the BLR), we also fit the spectrum with a power-law 
combined with absorption . In particular, absorption at E> 20 GeV/(l + z) is expected 



( ITavecchio &: Mazin 



20091 ) due to the Ly continuum, or at E>5 GeV/(l + z) due to the 



Hell recombination continuum ( jPoutanen fc Stern 1120101 ). We therefore fit the gamma-ray 
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Table 3. Gamma-ray spectral properties of GB6 J1239+0443, uncertainties at 68% C.L. 



Power Law 


Power Law + Double Absorber 


X z /DoF photon index 


X 1 /DoF photon index rjj rjfe 


1.1 2.14± 0.08 


1.6 2.13± 0.08 l-Oli'o °-0 9 



spectrum of GB6 J1239+0 443 with a power-law com bined with the gamma-gamma absorp 



tion model as proposed by 



Poutanen fc Stern 



(|201Q| ). Here, the absorption was fit with two 



parameters: (1) the optical depth for the Hi complex (r#) and (2) the optical depth for 
the Hell complex (r# e ). We fit the models to the data for gamma-ray energies below 20 
GeV, because the extragalactic background light is ex pected to absorb g amma-rays of E>20 
GeV, for sources at the redshift of GB6 J1239+0443 fjFinke et al. Ikoiol ). The results of the 
fit are reported in table [3j The fit with the gamma-gamma absorption components results 
in no absorption from the Hell complex, and weak absorption from the Hi complex. We 
performed the F-test on the two fits reported in table [3] to test the hypothesis of the need of 
the absorption components. The F-test gives a value of 0.15, and an associated probability 
of ~87%, hence the absorption component is not necessary to fit the spectrum, suggesting 
that in GB6 J 1239+0443 the blazar zone is located in the outer low- ionization region of the 
BLR or outside it. 



5.5 The spectral energy distribution for the high activity periods 

We constructed two spectral energy distributions referring to the high gamma-ray activity 
observed by AGILE and Fermi (campaigns A and B, respectively, as reported in Tabled]) 
with the data collected so far. The AGILE data are integrated for 4 days (the duration of 
the last observation during the campaign A, from 2008 January 4 13:35 to 2008 January 8 
11:06 UT). The Fermi-LAT data for period B are integrated within 4 days around the flux 
peak (2008 December 29 16:00 UT) to achieve acceptable statistics. These spectral energy 
distributions are shown in Figure [61 where we also included the archival data. 
The SED modeling has been performed in the framework of leptonic models. We assume 
that the emitting region is a spherical blob of radius Ruob, with bulk Lorentz factor T^fe, 
an electron population distribution proportional to 1+ ^j^_ sl+a2 (where 7 is the Lorentz 
factor of the electrons, ranging from 7 m j n to ^max)-, arid a randomly oriented magnetic field 
B filling the dissipation region. The observer line of sight and the jet direction form an angle 



®view We parametrize the external radiation fields as proposed by 



Ghisellini fc Tavecchio 
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Figure 6. In the two plots we report the SEDs for the two high gamma-ray activity periods detected with AGILE (red symbols, 
data from INTEGRAL/OMC, INTEGRAL/IBIS, AGILE-GRID) and Fermi-LAT (black symbols, data from VLBA, Kanata, 
Swift/XJVOT, Swift/XRT, Fermi-LAT) respectively. Green symbols are for data from Fermi-LAT integrated for one month 
period centered on the gamma-ray flare; light-cyan symbols are from the second Fermi-LAT catalog. We report optical data 
from SDSS (yellow) and UV data from GALEX (blue), as well as the Planck data (cyan), and the ROSAT detection (yellow) 
during the all-sky survey. Archival radio data from NED are shown in yellow. The disk model for the low activity period is 
plotted as a black dashed line; the model for jet emission in high activity period is reported as a continuous black line. In the 
top (bottom) plot we show the model assuming a blob dissipating at 0.2 pc (7 pc) from the central BH. 
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Figure 7. Same as previous figure, but in the jet emission model we assume a blob dissipating at ~5 pc, relaxing the condition 
Rfc!o6 = ^jR'di3si ancl assuming a blob radius of ~ 10 17 cm. 



(12009I ). where the key elements are the accretion disk luminosity and the distance (Rdiss) of 



the emitting blob from 



;he c entral BH. We consider a jet aperture 



Rblob 



Rdi. 



=0.1, according to 



Ghisellini fc Tavecchiol (120091 1 . In evaluating the power carried by protons, we assume one 
proton per emitting electron. We assumed that during the campaigns A and B, the accre- 
tion disk luminosity is the same as the one obtained from the SDSS data in 2001 (dashed 
curve in Fig. [6]). We found in section 5.4 that the Fermi-LAT gamma-ray spectrum (Fig. [6] 
green data, one month integration) does not show absorption from the BLR. Therefore, we 
tried to find solutions for the blob dissipating beyond the BLR, without taking into account 
the gamma-gamma absorption from the BLR. We obtained two possible solutions for the 
modeling. Model parameters are reported in Table HI where we use Rblr {Rtotus) the refer 
to the distance of the BLR region (of the dusty torus) from the BH; Iblr iJ-Torus) to the 
fraction of the disk emission that is reprocessed by the BLR in lines (by the dusty torus); 
taccr to the accretion efficiency; ^cooling to the electron Lorentz factor for which the electron 
energy halves in the blob crossing time (Rblob/ c). 

The first solution (model 1, top panel of Figure [6]) places the dissipation region (R^ss) just 
beyond the BLR (R^s ~0.2 pc from the SMBH), where the photon field from the dusty 
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torus is almost equal to the BLR seed photon field. At such a distance, the variability time- 
scale of a blob of radius Ruob ~7x 10 16 cm is of the order ~3 days, qualitatively in agreement 
(within a factor 2) with the 1-day binned gamma-ray light curve of the source. According 
to this modeling of the seed photon fields, the BLR photon contribution is expected to fade 
over time while the torus contribution remains constant. Both SSC and external Compton 
(EC) with seed photons from the torus contribute to the soft X-ray emission, the hard X-ray 
is dominated by EC emission with seed photons from the torus, and the GeV emission is 
dominated by EC with seed photons from the BLR. In the MeV up to GeV range the two 
EC contributions are almost equal. Further investigation of the variability is not possible, 
as no optical light-curve is available from which to trace the electron population from the 
synchrotron emission. Moreover the evolution of the gamma-ray photon index with time 
(reported in Fig. H]) is statistically poor. 

The second solution to the SED modeling (model 2) has been found by assuming Rrfj SS ~7 
pc from the SMBH. Such a model requires a Ruob = 2 x 10 18 cm, with a variability time-scale 
of the order of ~10 2 days. We do not find such a variability time-scale in the gamma-ray 
light curve, but we observe that the activity period for the source lasts at least 11 weeks. 
We can argue that the dissipating region in the jet has a radius of 2 x 10 18 cm. Gradients 



of the electron en ergy density in t 



re distance from the black ho 



(see for examples 



Giannios et al. 



2010 



Bromberg fc Levinson 



e, or disturbances in the jet 



20091 ) can be responsible for 



the shorter variability observed in gamma-rays. With this model, the X-ray emission is due 
to SSC; the gamma-ray emission is due to EC with seed photons from the torus. The direct 
contribution of accretion disk photons to the EC is negligible. 

In the model for a blob dissipating just beyond the BLR, the energy region below ~10 11 Hz 
is self absorbed. On the contrary, in the model of a distant dissipation region, there is no 
self-absorption, thus allowing the direct observation of radio emission from the blob. In this 
case, the emission predicted by the model is slightly higher than the flux reported in Fig. [6] 
(bottom panel), but we note that no radio observations were performed during the gamma- 
ray flares. The mismatch in the variability time-scale for model 2 with R^ ss ~7 pc, results 
from the assumption that Rbiob=j(jRdiss- Alternatively, we tried to model the SED assuming 
a blob of radius ~ 10 17 cm to satisfy the observed variability time-scale, and dissipating at 
~5 pc from the SMBH (model 3). Model parameters are reported in the right column of 
Table HI and the model is plotted in Fig. [7J 



A similar kind of solution has been suggested by iTavecchio et al. I (120111 ) to explain faster 
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Table 4. Model parameters for the fits of the spectral energy density, in the two assumptions of a blob dissipating just beyond 
the BLR (Rdiss ~ 0.2 pc) or far away from the SMBH (Rdiss ~ 7 pc). In the last column we report the model parameters for 
a blob dissipating far away from the SMBH (Rdiss ~ 5 P c )j but relaxing the condition Kuob= JqR<Hss> an d assuming a blob 
radius of ~ 10 17 cm. 





model 1 


model 2 


model 3 


R-diss(pc) 


0.22 


6.8 


4.8 


Blob radius (cm) 


6.7xl0 16 * 


2.1X10 18 * 


lxlO 17 


m S H (m©) 




5.3x10 s 




L d (erg/s) 




8.8xl0 45 




RflLfl (cm) 




3.0xl0 17 




R-Torus(cm) 




7.4xl0 18 




{blr 




0.1 




^torus 




0.3 





^■accr 0.1 



r bu i k 20 20 20 



angle of view (deg) 


2 


2 


2 




1 


1 


1 


'J max 


3.9xl0 3 


3.4xl0 4 


1.3xl0 4 


T&reafc 


0.95 xlO 3 


lxlO 3 


1x10 s 


density at -y break 
(cm -3 ) 


3.0xl0~ 2 


1.5xl0~ 4 


9.6xl0 -3 


81 


1.1 


0.5 


1.3 


S2 


3.1 


3.3 


2.5 


B (Gauss) 


6.1xl0 -1 


l.lxlO -2 


7.6xl0" 2 


^cooling 


60 


l.lxlO 4 


2.4xl0 3 


electron power 
(erg/s) 


2.2 xlO 45 


4.5 XlO 46 


l.lxlO 46 


magnetic power 
(erg/s) 


2.5xl0 45 


7.9 XlO 44 


8.6xl0 43 


proton power ** 
(erg/s) 


l.lxlO 47 


1.6xl0 47 


3.0xl0 47 


radiated power 
(erg/s) 


3.1 x 10 45 


2.5xl0 45 


2.1 x 10 45 



* R&iob=^)Rdiss in this model. 
** We assume one proton per emitting electron. 



variability (of the order of 10 min) for PKS 1222+216. This model can 
nection event, causing the acceleration of a small portion of the jet 



or to recollimation by the interaction with the external medium (jBromberg fc Levinson 



correspond to a recon 



Giannios et al. 



2010) 



2009|). 



According to this SED modeling, the X-ray emission is due to SSC, and the gamma-ray 
emission to EC with external photons from the torus. 



Flaring and low activity periods of GB6 J1239+0443 25 
6 DISCUSSION AND CONCLUSIONS 



We rarely have the opportunity to detect the disk emission in FSR Qs, which are gener- 
ally overwhelmed by synchrotron jet emission (see 



Pian et al. 



19991 . for example). How- 



ever, our study here suggests that we have detected accretion-disk-dominated emission in 
GB6 J1239+0443. Granted, we cannot fully exclude the possibility that the archival SDSS 
and GALEX observations we have reported could be interpreted as other emission mecha- 
nisms than thermal disk emission because we lack strictly simultaneous radio observations 
and extended radio light curves to corroborate the assumption of low jet emission. However 
the assumption of disk emission remains, in our opinion, the most likely explanation for the 
observations. 

The optical observations of GB6 J1239+0443 revealed an optical flux enhancement of a fac- 
tor 15-30 in 6 years, signifying a shift from accretion-disk to synchrotron, jet-dominated 
emission. The optical spectrum obtained in the period of faint optical emission allowed the 
classification of the source as a FSRQ in BZCAT. We made an estimate of the SMBH mass 
of GB6 J1239+0443 and the accretion rate from a period of low jet activity. With these esti- 
mates, we were able to study the December 2008 flare. Modeling the observed flat gamma-ray 
spectrum and SED also allowed for an investigation into the location of the blazar-zone of 
the object. 

As a final remark, it is worth stressing two major points. First, by definition, our estimate 



of Rdiss is model dependent. The 
the parametrization proposed by 



ocation of the dissipat i on reg ion was estimated assuming 



Ghisellini &: Tavecchiol ( 120091 ) for the BLR and the torus 



contribution to seed photons for the EC. According to this parametrization, it was possi- 
ble to derive Rdiss from the luminosity rat i o of sy nchrotron to EC emission. In fact, in the 
parametrization by iGhisellini fc Tavecchiol ( 120091 ) for a dissipation region outside the BLR, 
the seed photons for E C fade with distance from th e SMBH. We obtained two solutions. 



Referring to Figure 2 of 



Ghisellini &: Tavecchiol (120091 ). the ratio 



equals the ratio 



(U'blr+U'ir) 

of the optical to gamma-rays luminosity at the two values of the Rdiss (assuming knowledge 
of the magnetic energy density). For the flares of GB6 J1239+0443 reported in this article, 
one solution (model 2) places the blob at Rdiss ~7 pc, with only the dusty torus as the origin 
of seed photons. The other solution (model 1, with Rdiss ~0.2 pc) corresponds to just outside 
the BLR, where the contribution of seed photons from both the dusty torus and the BLR are 
relevant. The magnetic field is constrained in the SED modeling by the cutoff of synchrotron 
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emission in the UV (due to the last and most energetic electrons), and by the corresponding 
cutoff of EC which we cannot derive directly from data (that give non-constraining upper 
limits at E>20 GeV, see Fig. [6]). Assuming the Thomson regime, and with only one external 
photon field contributing to the EC, the ratio f cu toff between the synchrotron cutoff energy 
and the EC cutoff energy is proportional to r& ^ u — — , where < v seed > is the typical seed 
photon field energy. Hence if we have constraining data at the highest energy, and with a 
specific geometry in the model, we can constrain ~B/Tb u ik- However, the geometry in the 
model also constrains U' BLR + U' IR , hence in the ideal case we can obtain B/r^ and Rdiss 
from the SED modeling. The ratio of synchrotron to SSC luminosity further constrains the 
model parameters Rbiob,^buik, B, allowing one to remove the degeneracy between T bu i k and 
B. In reality, we can obtain only upper limits of f cu toff because the data at higher energies 
are non-constraining. So we have only upper limits on B/T^ik for each model. As a conse- 
quence, the U' BLR + U' IR could be lower than in our parametrization (we have to maintain 
u' 

the ratio 777 — ttp~ at the desired value). This implies that Rdiss could be higher than our 

BLR I R 

evaluations. 

The second point is that the photon field intensity is proportional to the accretion disk lu- 
minosity, and the BLR and torus location is proportional to VXd- in all our estimations, we 
assume that the disk luminosity is almost steady over time, e.g., in the low state observed 
during the Sloan Survey in 2002, during the GALEX observation in 2007, and during the 
gamma-ray flares observed by AGILE at the beginning of 2008 and by Fermi-LAT at the 
end of 2008. This assumption could be false. In this case, we observe that t h e par ametriza- 
tion of U' BLR and U' IR reported in equation (20) by iGhisellini &: Tavecchiol ( 120091 ) remains 
unchanged while varying L^, provided that we scale the solution for Rdi SS with \[L~d- There- 
fore, variations of the disk luminosity in time, and/or systematic errors in the evaluation of 
disk luminosity from our SDSS+GALEX data (possibly biased by jet emission) only slightly 
affect our estimation of Rdiss- 

The starting point of our modeling is that the emission region is far from the SMBH (at 
parsec scale), and we motivate this choice with the flat gamma-ray spectrum up to energies 
of 15 G eV. For a different app roach and results for other blazars, we refer readers to the 



work of 



Tavecchio et al. 



( 120101 ). They performed a detailed study of the localization of the 
emission region for bright blazars making use of the variability timescale for objects showing 
a spectral cut-off at (10-20)/(l+z) GeV. They obtained that the variability timescale and 
spectra are in agreement with a dissipation region inside the BLR. 
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We note that, contrary to the model we used, some authors (e.g., 



Giommi et al. 



2011 



, for 



example) model the SEDs of both FSRQ and BL Lacs with pure Synchrotron + Synchrotron 
Self Compton components only. 

Model 1 (Rdiss ~ 0.2 pc), gives a variability time scale of the order of 3 days, in agree- 
ment with the flare duration estimated from the 1-day binned gamma-ray light curve. This 
model, however, does not properly reproduce the flat gamma-ray spectrum. In particular, in 
order to reproduce the observed flux at energies > 10 GeV, it overestimates the spectrum 
for lower energies. On the contrary, model 2 (Rdiss ~7 pc) reproduces the flat-gamma-ray 
spectrum, but it is not clear whether the predicted variability of the order of ~100 days 
can be associated with the duration of the gamma-ray activity period estimated from the 
1-week/ 1-month binned gamma-ray light curves. The third model has been built by relaxing 



the relation Rbi o b=0.lRdi S s in order to preserve the variability tim e-scale estimated 
1-day gamma-ray light curve (we follow the solution proposed by 



Tavecchio et al. 



rom the 
2011 for 



PKS 1222+216), and it still reproduces the gamma-ray spectrum. Interestingly, the size of 
the emitting region for PKS 1222+216 (R Uo b ~ 5 ■ 10 14 cm) and for GB6 J1239+0443 differ 
significantly. 

For the third m odel, we obtain R h i nh =0. 067 R diss, in agreement within a factor of two with 
the prediction of iBromberg fc Levinson I (120091 ). that gives Ru b — 10 _2 ' 5 -R* SS , for the case of 
efficient conversion of bulk luminosity in radiation in the strong focusing scenario. With the 
same assumptions, IBromberg &; Levinson I (120091 ) assume that the location of the emitting 
region is at Rdiss ~2.5 (Lj e i/10 46 erg s~ 1 )(Rblr/0.1 pc)~ x pc from the SMBH, where Lj et is 
the jet power. If we invert this relation, and we make use of our result (Rdiss = 4.8 pc), we 
obtain Lj et ~ 3.5 • 10 46 'ergs -1 . We must assume that the proton-to-emitting electron ratio 
is of the order of 0.1 in order to reproduce such a power (in the evaluation of proton power 
reported i n Table IH we assumed one p roton per emitting electron, instead). We note, how- 



ever, that iNalewajko fc Sikora I (120091 ) evaluated that efficient radiative conversion could be 
assumed if the product of the bulk Lorentz factor by the opening angle is >3, and according 
to our third model this product is 2. 
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